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Abstract 

We have searched for long periodicity in ten X-ray 
selected Blasers with RXTE ASM 4.6 years data. 
We found about 10-100 day possible periodicities 
for three TeV gamma ray Blasers, Mkn 421, Mkn 
501 and PKS 2155-304. 



using the RXTE ASM archive data with a Lomb 
method. 

2 Observation 

We used RXTE ASM data from January 1996 
(MJD 50087) to August 2000 (MJD 51780), (4.6 
years). We collected Xray selected Blasers cata- 
loged by Xie et al.(1993) from ASM archives, as Ta- 
ble 1. We obtained a series of 90 sec integrated data. 
We show the lightcurve of three interesting Blasers 
which show possible periodicities in our analysis, in 
figure 1. 



1 Introduction 

Begelman et al.(1980) suggested that multiple 
Blackholes may exist at the center of AGN from the 
existence of merging galaxies. If there is a multiple 
Blackhole system, we may see an orbital period by 
a secondary Blackhole crossing the accretion disk 
of the primary or a precession of a jet. For exam- 
ple, a 12 years periodicity is confirmed in optical 
data with 100 years for a radio selected Blaser OJ 
287(Sillanpaa et al. 1988). Also, a periodicity of 
23 days for X-ray selected Blaser Mkn 501 was de- 
tected with the TeV gamma ray cherenkov detec- 
tor Telescope array (Hayashida et al. 1996) and 
HEGRA (eg. Aharonian et al. 1999a) and with the 
Xray detector RXTE All Sky Monitor (ASM;Levine 
et al. 1996) during a TeV flare in 1997 (Kranich et 
al. 1999; Nishikawa et al.l999). 

There is an observation of radio jet bending in 
Xray selected Blasers Mkn 421 (Finer et al.), Mkn 
501 (Giovannini et al. 2000) with VSOP, which 
may suggest precession or orbital motion of the jet. 
There also may be a long periodicity from a thermal 
instability in an accretion disk (Honma et al. 1991; 
Abramowicz et al. 1995). If we scale the periodicity 
of 10 sec of Blackhole binary GRS 1915-^105 (eg. 
Morgan et al. 1997) to AGN of a mass with lO'^ 
Mq, a periodicity of 100 day is expected. 

Here, we concentrate on Blasers which may have 
a periodicity from a geometrical origin. We aim to 
confirm a periodicity of 23 days for Mkn 501 and 
search for longer periodicities than 1 day for Blasers 



Table 1 : Target list of Xray selected Blasers 

Mkn 501 FKS 0548-322 Mkn 421 

Mkn 180 H 1101-232 BL 1ZW187 

BL 1426-H427 BL 12194-305 BL 0323+022 

FKS 2155-304 



3 Analysis 

In order to examine the periodicity of X-ray intensi- 
ties from each source, we use a Lomb method (Lomb 
1976; Scrargle 1982), which is suitable for the anal- 
ysis of an unevenly time coveraged data set. 

3.1 Data selection 

We used three detector data of ASM in order to in- 
crease the number of data points. Since we searched 
for a periodicity longer than 1 day, we did not carry 
out the baricentric correction. 

We selected data points with a significance: the 
ratio of the rate to error, S — rate/error, with lev- 
els of > 0,1,2,3,4 in order to control the data qual- 
ity. As well as raw data (no binned data), we also 
binned a 90 sec integrated data in 1 day, 5 day, 10 
day intervals to solve poor statistics. In binning, 
we used a weighted mean and obtained errors by 
summing in quadrature. The power in the Lomb 
method is linear in the number of data points if 
a periodicity is real. As a result, by binning or 
higher S selection, we obtained a fewer number of 



(lata points and a lower significanc;e for a pcriodic;- 
ity. Therefore, we checked for the periodicity under 
all conditions mentioned above. Even after binning, 
some data points arc still statistically poor because 
the original data point distribution is not uniform in 
time. Therefore, we selected data points consisting 
of > 20 raw data points. We searched for periodicity 
with a data set of 1700 days. For Mkn 501, we also 
search periodicity with a data time limited to MJD 
50300-50900, which corresponds to time duration of 
a long TeV flare. 

3.2 Result 

We calculated a power spectra with a number of 
frequencies of lOOxA^/2 (iV:number of data). We 
limited the search interval for periodicity between 
1 day and Tohs/10. Here, Tots is the total observa- 
tion time. Wo calculated a chance probability with 
Prob(> 2) = 1 — (1 — e~^)^ . Here 2; is a power. We 
detected possible periodicities with a chance prob- 
abihty less than 10~^ for Mkn 501, Mkn 421, and 
PKS 2155-304. 

3.3 Consideration of a gradual in- 
crease 

When we applied the fourier transformation or more 
sophisticated Lomb method to the limited time in- 
terval data, we always saw spurious large powers 
around the test period equal to the observation 
time. This is due to the low frequency fourier com- 
ponent of Xray data. We have confirmed this ef- 
fect by simulation. For example, we have gener- 
ated a linearly growing lightcurvc in the observa- 
tion time, rate=a*t -|- b with error sampled from 
real data, then we analyzed this data with the 
Lomb method and found a similar spurious peri- 
odicity around P ~ Tobs- Therefore, to exclude 
such effects, we removed the low frequency com- 
ponent as follows. We fitted the lightcurvc with 
a polynomial function with order n=l,2,3,4,5, then 
we employed the best model, which gave a mini- 
mum value of the x^/d.o.f, as a global lightcurvc. 
The refined lightcurvc obtained by subtracting the 
modeled function (global lightcurvc) from the orig- 
inal lightcurvc was analyzed. Furthermore, we lim- 



ited the periodicity search range between 1 day and 
To6s/10 to reject any spurious periods. As seen in 
figure 1, the lightcurvc of Mkn 421 shows double 
peaks with separation of 900 days. Even after sub- 
tracting model function, we could not remove the 
interference between these two peaks and we ob- 
tained very large power at 900 days period. There- 
fore we only used the Mkn 421 data between MJD 
50088 (beginning of observation) and MJD 51500 
(just before the second peak). 

We show the results in table 2 and figure 2. 
The many periodicities which have probability with 
more than 10~^ in figure 3 are considered to be sta- 
tistical fluctuations. 

4 Discussion 

For Blasers, there is no origin for a long periodicity 
in a Blackhole system. However, if multiple Black- 
holes exist in the center of AGN, we can see a pre- 
cession or an orbital motion by a secondary Black- 
hole crossing an accretion disks of the primary or 
by orbital (j) component. Hence, we consider three 
Blasers with significant periodicities as candidate 
for multiple Blackhole systems. We discuss param- 
eters of multiple Blackhole systems with the infor- 
mation of the obtained periodicities and discuss a 
relation between a periodicity in Xrays and TeV 
gamma rays. 

4.1 Geometry with orbital period 

We consider a two Blackhole system for simplic- 
ity. Assuming a detected periodicity as an orbital 
period, we estimate a separation distance of two 
Blackholes as shown in table 3 using the formula 
for Kepler motion, 

QM1+M2 _ (22:^^ 

Here, Ml, M2 is each Blackholes mass, P is the or- 
bital period, r is the separation radius. Even with 
a current high resolution instrument such as VSOP 
( angular resolution 0.1x10"^") or a next genera- 
tion VERA ( 10x10"^"), it is not possible to resolve 
such a two Blackhole system. We hope for a confir- 
mation of multiple blackholes by a high resolution 



Table 2: Periodicity for all selections. The periodicities which have more than 10 cycles and a chance 
probability less than 10~^ is shown. Low frequency component has been removed. 



Target 


binning(day) 


S 


Period (day) 


Chance probability 


Number of data 


Mkn 421 


no 





62.1 


8.7x10"^ 


18912 








107.4 


8.2x10-5 




Mkn 501 


no 





47.4 


5.8xl0-« 


27678 


Mkn 501(MJD50300-50900) 


no 





23.6 


2.2x10"*^ 


9778 




no 


1 


23.6 


3.9x10-4 


5157 


PKS 2155-304 


no 





142.9 


5.1x10-^ 


17897 




no 


1 


143.3 


4.4x10-'^ 


7166 



Table 3: Calculated geometry, assuming the de- 
tected periodicities caused by an orbital period. 
Ml + M2 = 2x10^Mq is assumed. (i?o = 65 
km/(s-Mpc)). 



Target 


period 


separation 


angular 




(day) 


radius (pc) 


distance 


Mkn 421 


62.1 


9.0x10-4 


1.3xl0-«" 


Mkn 501 


23.6 


4.6x10-4 


6.0x10-^" 


PKS 2155-304 


143.0 


1.5x10-3 


6.1x10-^" 



instrument sometime in the future. 



4.2 Constraint on mass ratio 

A precession period for a differentially rotating fluid 

disk is given by Papaloizou and Terqucm (1995) and 
the precession period Pprec is related to the orbital 



period P by Larwood (1998) as, 

Pprec 7^^\1+IJ.) V 



V 3/2 



cos6. 



Here, Pproc is the primary precession period and /i 
is Mg/Mp, Mp is the primary blackhole mass, Mg 
is the secondary one, r is the separation distance 
between the two Blackholcs, i?disk is the radius of 
the accretion disk, and 6 is the orbital inclination 
angle with respect to the disc. 
Now, we limit the mass ratio between two black- 
holes from the detected periodicities, using the 
above formula. We assume the observed periodic- 



ities as precession periods Pp., 



The two Black- 



hole system with a short orbital period can not 
be alive for a long time because of gravitational 
wave emission. When we assume a typical Black- 
hole separation radius r/i?disk '^100, the condition 
of ~ 10^ gives an acceptable orbital period of 
:0.14 cosi5 Pprec ~ 20 days for Pprec=140 days. If 
we assume the condition with smaller /i, for exam- 
ple /x=l, we obtained unacceptable orbital period 
of 3x10-4 cosSPprec- We conclude the mass ratio 
should be very large, typically ^ > 10'''. 
On the otherhand, if we assume observed periodic- 
ties as orbital periods, there are two possible mod- 
els. One is that we can sec the orbital periods by 
a secondary crossing the accretion disk of the pri- 
mary or two accretion disk causing a tidal distur- 
bance. This needs a large mass difference for two 
Blackholes (Sillanpaa et al. 1988;Letho & Valto- 
nen 1996). The other is that we can see the orbital 
period by the orbital (f> motion, which needs two or- 
ders mass difference for 23 day periodicity of Mkn 
501(Rieger & Mannheim 2000). 

4.3 Relation with TeV 

It is very interesting that the three Blasers 
with siginificant periodicities are all TeV gamma 
sources (Mkn 421;eg.Aharonian et al. 1999b, Mkn 
501;eg. Aharonian et al. 1999a, PKS 2155- 
304;P.M.Chadwick et al. 1999). If these period- 
icities are due to a geometrical effect, we expect 
a similar periodicity in TeV gamma ray intensi- 
ties. We searched for such a periodicity using a 
Lomb method with HEGRA TeV published data, 
for Mkn 421 during 1997-1998 and for Mkn 501 



during 1997-1999 (Loronz 1999; Aharonian 1999b; 
Aharonoan 2001;Kranich 2001). We confirmed a 23 
day periodicity for Mkn 501 during the 1997 flare 
time found by Kranich ct al.(1999). However, dur- 
ing 1998-1999, 23 day periodicity was not found for 
Mkn 501. A 62 day periodicty was not found for 
Mkn 421. Probably, the sensitivity of TeV gamma 
ray Cherenkov detectors may not be enough to de- 
tect periodicities except around the large flare of 
Mkn 501 in 1997. 

The good coincidence between periodic Blasers and 
TeV gamma ray Blasers is remarkable. We can es- 
timate the chance probability of this coincidence at 
less than 10~^. There may exist a relation between 
multiple Blackholes and an electron acceleration to 
ultra high energy. 
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Figure 1: Lightcurve of three interesting Blasers in our analysis. A possible periodicity with chance 
probability less than 10~^ is detected. Vertical axis is the count rate for each detector of ASM , horizon 
axis is the observation time in unit of MJD. Data is binned in 50 day increments and selected with data 
points >20 for clarity. 
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Figure 2: Power spectra for Mkn 421, Mkn 501, Mkn 501 during MJD 50300-50900 (around the TeV 
flare time), PKS 2155-304. Low frequency components have been removed. This data has been obtained 



